function y0=Neville_eval(~,Y,~)
%X，Y是已知的插值点坐标点%x0是插值点
% y0是多项式在x0处的值m=length(X);
P=zeros(m,l);
P1=zeros(m,1);
P=Y;
for i=1:m
P1=P;
k=1;
forj-i+l:m;
k=k+1;
P(j)=P1(j-1)+(P1j)-P1(j-1)*(x0-X(k-1))/...X(j)-X(k 1));
  end
if abs(P(m)-P(m-1))<10-6
y0=P(m);
return;
end
end
y0=P(m);